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The matrix model version of AGT conjecture and CIV-DV prepotential 

A. Morozov 1 and Sh.Shakirov 2 
ABSTRACT 

Recently exact formulas were provided for partition function of conformal (multi-Penner) /3-ensemble in the Dijkgraaf-Vafa 
phase, which, if interpreted as Dotsenko-Fateev correlator of screenings and analytically continued in the number of screening 
insertions, represents generic Virasoro conformal blocks. Actually these formulas describe the lowest terms of the g a -expansion, 
where q a parameterize the shape of the Penner potential, and are exact in the filling numbers N a . At the same time, the older 
theory of CIV-DV prepotential, straightforwardly extended to arbitrary f} and to non-polynomial potentials, provides an alternative 
expansion: in powers of N a and exact in q a . We check that the two expansions coincide in the overlapping region, i.e. for the lowest 
terms of expansions in both q a and N a . This coincidence is somewhat non-trivial, since the two methods use different integration 
contours: integrals in one case are of the B-function (Euler-Selberg) type, while in the other case they are Gaussian integrals. 
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1 Introduction 

This paper is a part of extended study of AGT conjecture [1] - [6], which unifies a number of important subjects 
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[8] - [23], that were in the focus of theoretical investigations during the last decade. This particular paper 
is devoted to the recently formulated explicit relation [4, 5, 6] between the integrated free-field correlators in 
Virasoro conformal field theory (CFT) on the Riemann sphere 



Zdf = Jd Zl ... J dz N /jj : e 5 «^«) : J] 



e b<i>(zi) 



C± Cn 



= no** - x «? &a&b J dzi... j dzvUizj - Zl ) 2t > 2 y[( Zi - x a ) 2ba ° (i) 

a<b Ct Cn i<j i,a 

and conformal blocks - holomorphic and anti-holomorphic parts of generic (not necessarily free-field) CFT 
correlators. Actually, this relation was long expected: since the old papers [7], there was a belief that any 
spherical Virasoro conformal block [8] can be obtained by appropriate choice of integration contours C\ , . . . , Cn 
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in the free- field Dotsenko-Fateev integral Zdf- However, precise formulation (including the choice of contours) 
remained a mystery for some time and was made explicit only recently [5]. 

The aim of the present paper is to investigate the partition function Zdf by the methods of matrix models. 
As noticed long ago [22], integrated free-field correlators like (1) can be viewed as natural objects in matrix 
model theory, which studies ensembles of " eigenvalue" variables z\, . . . ,Zn with a certain power of Vandcrmondc 
determinant in their integration measures 



N 

:.. , F ~ / dzt . . . I dz N Y[( Zj - z t ) 2 P exp ( - \, I3 = b 2 (2) 

i<j V i=l / 



Ci Cn 

where in the case of Zdf potential V(z) is of the multi-Penner shape: 

r+l 

V(z) = -^2a a log(z - x a ), a a = 2ba a (3) 

o=l 

Models with such a /3-dependent eigenvalue measure are known in literature as /3-ensemblcs [21]. For generic 
values of j3 they are not, strictly speaking, matrix models, since they cannot be induced by integration over 
any kind of matrices. However, it appears that most of conventional matrix-model methods and technical tools 
(such as Ward identities and genus expansion) are perfectly applicable for arbitrary values of /3. 

Because of this, Dotsenko-Fateev partition functions can be viewed as (and treated essentially by the same 
methods as) non-Gaussian eigenvalue matrix models [13, 15], only taking into account properly the /3-deformed 
integration measure. In this paper we consider a pair of alternative approaches to such non-Gaussian integrals: 
1) based on the quasiclassical expansion and 2) based on exact contour integration. Our ultimate aim is to 
compare the results of these two approaches with each other and elucidate all possible subtleties. 

1) The first method we use in this paper is quasiclassical approach to matrix models: namely, the theory 
of Dijkgraaf-Vafa (DV) phases [13, 15]. This approach exploits the fact, that non-Gaussian integrals 



J dz 1 . . . J dz N Y[(zj - Zi) 2f} exp (^-^V(z 1 , . . .,z N )j 



(4) 



in the limit h — > arc dominated by fluctuations of integration variables near one of the saddle points 
(Zi, . . . , Zn) - solutions of the saddle point equation (or, what is the same, equation of motion): 

Since for the eigenvalue matrix models the full multivariate potential has a form of a sum V{z\, . . . , zn) = 
V(zi) + . . . + V(zjy), each of components Z^ has to be a critical point of the single-variable potential V(z). 
Accordingly, the N components of the saddle point Z\, . . . , Zpj get divided into r groups: 

(Z 1 ,...,Z N ) = (//i,...,/ii,/Z2,...,/i2,...,Mr,---,Mr) (6) 
s v ' * v ' * v ' 

ATi N 2 N r 

where fii, . . . , fi r denote the critical points of V{z) (i.e, solutions of V'(fi a ) = 0) and r denotes their total 
number. The parameters A^ satisfy Ni + . . . + N r = N and label the choice of physical phase for the multi- 
Penner /3-ensemble. After the choice of particular set {N\, . . . , N r }, one needs to integrate over all fluctiations 

(z!,...,z N ) = (Zx,...,Z N ) + (yi,...,yjv) (7) 

treating the corresponding integrals over fluctuation variables y\, . . . , t/jv as Gaussian correlators. This method 
is quite general and applicable to generic potentials; applying this method to particular non-Gaussian integral 
(1), one obtains an answer for (1) in terms of power series in the small parameter h. 

2) The second method is to specify concrete integration contours for the non-Gaussian integral and 
integrate over these contours directly. Sometimes called theory of integral discriminants [23], this approach 
is rarely used in the literature: since the methods of exact integration for non-Gaussian potentials are under- 
developed, their application in practice is feasible only for simple enough potentials V(z). Fortunately, the 
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potential V(z) we are interested in - the multi-Penner potential - has a very special form: it is a logarithm of a 
polynomial (for all natural parameters a a ). Consequently, the exponent exp ( — V{z)) is merely a polynomial, 
so that exact non-Gaussian integration for the partition function (1) is feasible in practice. 

At the first glance, the choice of integration contours C\, . . . , Cn seems unrestricted, since it is quite easy 
to integrate a polynomial along any contour in the complex plane. However, as usual in contour integration, 
homotopic choices of contours actually give one and the same answer for the integral. Because of this, the set of 
essentially inequivalent choices consists of all segments j ah — [x a ,x ], which connect the roots of the integrand. 
Moreover, the obvious relations -f ab + 7 fcc = "f ac leave only r inequivalent contours j a — \x\, x a ] at our disposal, 
so that the N contours C\, . . . , Cn get divided into r groups: 

(Cl,...,C N ) = (72,...,72,73,---,73,---,7r+l,---,7r+l) (8) 

V v ' V v ' * v ' 

JVi JV 2 N r 

Note, that for any natural N ai a a and j3 the integrand is a polynomial, hence the integral can be evaluated 
exactly. Note also, that (just like in the first approach) additional parameters N a appear naturally. They satisfy 
N\ + . . . + N r = N and label the choice of phase for the multi-Penner /3-ensemble. One concludes that this 
appearance is inavoidable, reflects a natural property of non-Gaussian integrals (they are not fully defined by 
the action alone) and should be reproduced by any other consistent approach to non-Gaussian ensembles. 

Applying this method to the non-Gaussian integral (1), one obtains an answer for (1) in terms of power 
series in the small parameters q a (related to x a by a simple law, see [5]). These series were first evaluated in 
[5, 6] and found to coincide with the standard conformal blocks in Virasoro CFT, thus providing an explicit 
and convincing check of the "matrix-model version" of the AGT conjecture. The observation of [5] has lately 
led to a number of developments [6], aimed at proving the original AGT conjecture by matrix-model methods. 
This paper can be considered as yet another contribution to this direction of research. 

1) vs. 2) In the present paper, we adress the issue of the relation between two different methods 1) and 

2) , and two corresponding expansions - in powers of h (the quasiclassical DV expansion) and in powers of q 
(the CFT expansion). Both methods can be applied to calculate the partition function of the multi-Penner 
/3-ensemble Zdf- Despite they are designed to calculate one and the same quantity - partition function of the 
non-Gaussian /3-ensemble - it is not a priori obvious that their results are consistent with each other. 

Actually, the reason for doubts is severe difference in the choice of integration contours. The first method 
(selection of contours and exact non-Gaussian integration over these contours) naturally incorporates the choice 
7a = [xi, x a ] of integration contours for the variables zi,..., zn- Such contours are not quite typical for matrix- 
model theory: usually, one integrates either along closed contours or along open contours that go to infinity, like 
the real axis in the Hcrmitian matrix model, the imaginary axis for the anti-Hermitian model, etc. As explained 
in section 4, precisely such "typical matrix-model" contours are used in the second method for integration over 
the fluctuation variables yi, . . . , xjn- It is not a priori clear, why such different treatments of a non-Gaussian 
integral would give one and the same answer. 

Therefore, consistency between the two methods is not quite trivial and deserves to be checked. It is a goal 
of present paper to perform such a check. To simplify our presentation, we consider only the simplest case of 
the 3-Penner potential (which corresponds on the CFT side to the 4-point conformal block on a sphere): 

V(z) = -a± log(z) - a 2 \og(z - q) - a 3 log (2 - 1) (9) 

The first method provides a series expansion in a single small parameter q, described in [5, 6]. The second 
method provides a series expansion in h, derived in this paper. To compare both methods, we study the double 
(h, g)-expansion from both the contour integration and quasiclassics sides. We find complete agreement. 

This paper is organized as follows. We begin in section 2 with a simple (perhaps even oversimplified) 
example of comparison between 1) the contour-integration and 2) the quasiclassical DV methods. This example 
illustrates and helps to clarify our claims. Then, in section 3, we describe the DV method for arbitrary non- 
polynomial potentials in full generality. The particular case of logarithmic 3-Penner potential is considered in 
section 4: using the DV method, we obtain /i-expansion of the 3-Penner partition function and derive its double 
(h, q) expansion. Finally, in section 5, we consider the g-expansion of the same 3-Penner partition function and 
derive the double (h,q) expansion. A priori different, the two double (h,q) expansions are found to coincide. 
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The Appendix of the present paper is devoted to another important special case of the general considerations 
of section 3: namely, the case of polynomial potentials 



r+l 

V{z) = Y,T a z a (10) 
0=1 

Models with such potentials recieved a considerable amount of attention [13] some time ago, because at (3 = 1 
their free energy F — log Z makes appearance in supersymmetric 4d gauge theories with spontaneously broken 
gauge symmetry. The quasiclassical expansion, which we use most often in the present paper, is actually not 
the most convenient for comparison with 4d gauge theories: more convenient is the so-called genus expansion, 
that is, the same quasiclassical expansion expressed in terms of different variables S a — hN a . In terms of 4d 
gauge theory, variables S a posess an interpretation of vaccuum condensates of gauge fields, while variables N a 
do not literally correspond to any gauge theory quantities. In terms of variables S a , the free energy takes form 

= 1: F = logZ=^^F (S 1 ,...,S r ) + h 2 F 1 (S u ...,S r ) + ... ) (11) 

where the components F p (Si, . . . , S r ) are called genus-p prepotentials, and the leading genus zero component 
Fo(S\ , . . . , S r ) is widely known as CIV-DV prcpotential. Using the general formulas of section 3, we calculate 
the CIV-DV prcpotential (for arbitrary degree r = deg V(z) of the potential) and generalize it in two directions: 
to higher genera and to arbitrary ft. For j3 7^ 1, we observe appearance of contributions of odd order 

0^1: F = logZ= p( F (S 1 ,...,S r ) + hF 1/2 (S 1 ,...,S r ) + h 2 F 1 (S 1 ,...,S r ) + ... ) (12) 

which correspond to non-integer genera in the genus expansion. We derive explicit expressions for the generalized 
/3-deformed CIV-DV prepotentials of genera 0, 1/2, 1 and 3/2, and study some of their simplest properties. 
Despite these expressions have little in common with the main topic of present paper, they are interesting by 
theirselves and will be probably useful in future research. We include them for reference reasons. 

2 A toy model: B-function 

To illustrate the difference between the choices of integration contours in methods 1) and 2), let us now consider 
an oversimplified example: the B-function integral 



B(u,v) = J z u (l - zfdz (13) 

This integral, actually, corresponds to the particular case a.\ — u, a 2 — 0, a 3 = v and N = 1 of the 3-Penner 
integral. From the point of view of the method 1), as explained above, we need to choose essentially independent 
integration contours. We can choose one such contour - connecting the points and 1. With this choice of 
contour, the integral becomes correctly defined and convergent, giving an exact answer 



Bi(«,«) = J z u {l-zYdz=^^ 



r(tt + i)r(v + i) 

+ V + 1) 



(14) 



As already mentioned above, this contour is not typical for matrix models; a typical matrix-model contour 
would be (—00, 00). From the point of view of the method 2), we need to find the critical points of 

V(z) — — ulogz — wlog(l — z) (15) 

i.e, solutions of the equation of motion 

V{z) = -±-^ = (16) 
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It is easy to see that the only solution (the only critical point) is u/(u+v). To describe quasiclassical fluctuations 
of the single variable z around this critical point, we need to make a change of variable 

z = ^- + y (17) 

u + v 

where y describes the fluctuation. The DV procedure, just like any other quasiclassical method, prescribes to 
integrate the fluctuation from — oo to +00. Accordingly, the integral takes form (we omit ft for simplicity) 

00 

„ . , u u v v f ( (u + v) 3 2 (u + v) A (u-v) o \ , 

B > ^ v ^ - J cxp {-^r y - 3^ - y + -) dy (18) 

—00 

The point is that a priori Bi(u,v) and B 2 {u,v) are different integrals, since the contours of integration differ. 
However, they do coincide as expansions in u,v >> 1 (equivalently, as expansions in ft — > 0). One can check 
this directly by evaluating B2 (u, v) as a Gaussian integral 



_ , » u u v v / 2uv 1 — / v 2 — ' km -.- (/ \ 



(u + v) u+v V (u + v) 3 \ 12uv(u + v) 

what reproduces the well-known Stirling expansion of Bi(u,v) at u, v » 1. In fact, it is hard to expect a 
different result in this toy model: we consider just a quasiclassical expansion of a single- variable integral. What 
we check in present paper, is that the coincidence holds further for N > 1 and a 2 7^ 0. 

We now mention another way to derive the relation between Gaussian and B-function (hypergeometric) 
integrals, which may be more straightforward but somewhat delicate as soon as one touches the issues of 
integration contours and analytical continuation. From the definition of the exponent we have 



/oo roo / 2 \ Ivl poo 

z 2p e- z ' 2 dz= lim / z 2 p ) dz= lim M" +1 / 2 ^"^(l - X ) M dx 

If the integral at the r.h.s. was between to 1, then the r.h.s. would be 

lim M^(M + !)!>+ 1/2) = = (2^-1)11^ 

m^oo r(M + p + 3/2) ^ ' ' 2Pp\ 



(20) 



(21) 



what is the correct answer for the original Gaussian integral. If integral at the r.h.s. of (20) was indeed taken 
from to 00, it would diverge, but one could split it into two, from to 1 and from 1 to 00. The latter one 
could then be analytically continued with the help of modular transformation x — > 1/x: 



/OO f-\ 
x p -V 2 (l - x) M dx -+ (-) M j x- m -p-^ 2 (1 - x) M dx = 



r(l/2-p) 

(-) M sin(p + 1/2)tt T(M + l)T(p + 1/2) _ ' : 



(-) M T(-M-p-l/2)T( M+l) = 



sin(M + p + 3/2)7r T(M+p + 3/2) J 



x p -^ 2 (l-x) M dx (23) 



so that adding the two pieces one would get a vanishing answer. Of course, the right answer is to declare 
that the hypergeometric integral at the r.h.s. of (20) is a B-function, i.e. goes from to 1 only — and 
this is exactly what we observe in this paper in more general situation. In the simple context of (20) 
one can easily justify this prescription. Indeed, let original Gaussian integral be taken between — K and K with 
K ^> 1 but finite: this is a good approximation, of course. Then the x-integral in (20) will go between and 
J K/M, and this upper limit tends to zero if M — > 00 with K fixed. Thus there is definitely no need to include 
the region 1 < x < 00 into the integration domain. On another side it can be continued from \J K/M <C 1 to 
1, because for large M only the narrow region near x = actually contributes. 
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3 /i-expansion for the generic case 

For arbitrary potential V(z), one can define and study the non-Gaussian partition function 



Z= fdzt... ( dz N Y[( Zj - z^ exp l~ Yl ) ( 24 ) 

J J i<j \ i=l J 

Such models have been studied for quite a long time: in the case, when V(z) is a polynomial, their free energy 
F = \ogZ gives rise to the CIV-DV prepotential [13]. Throughout this section, we do not assume that V(z) 
is polynomial; we assume only that V(z) admits an infinite Taylor series expansion near each of its critical 
points. In particular, all polynomials and all logarithmic (multi-Penner) potentials do fall into this class, while 
potentials of fractional degree (like (z — ^) 3/ " 2 ) do not. 

It is most important to note that, since the integral is non-Gaussian, it is not completely defined by eq. 
(24): in addition, it depends on the choice of integration contours. There are at least two different ways to 
remove this ambiguity: either to select a particular contour (actually, one can associate an independent contour 
d with every variable Zi) or to study the integral quasiclassically around some particular saddle point. This 
section is devoted to the second - quasiclassical - option. 

The saddle point analysis of the integral (24) is quite simple: when ft — > 0, variables z, become distributed 
between the critical points - solutions of the equation V'(z) = 0. For arbitrary V(z), this equation can have 
infinitely many solutions. Two classes of potentials, often encountered in practice, posess a finite number r of 
critical points [i\ , . . . , \x T : polynomial potentials 

r+1 

V{z) = Y. T * z<1 ( 25 ) 

a=l 

and logarithmic (multi-Penner) potentials 

r+1 

V(z) = -J2a a log(z-x a ) (26) 

a=l 

It is easy to see that different saddle points of the non-Gaussian integral are labeled by different partitions 
N = Ni + . . . + N r . Indeed, the only possibility modulo permutations is that the first Ni eigenvalues tend to 
the critical point fi\, then next iV 2 eigenvalues tend to fj,2, and so on: 



(Z lt ...,Z N ) = (^i,...,/ii,...,/x r ,...,/x r ) (27) 

as illustrated on Fig.l. The numbers Ni are called filling numbers. To evaluate the integral quasiclassically, we 
need to study the integrand in the vicinity of the saddle point: 

(zi,...,z N ) = (^i,...,/ii,...,/i r ,...,/i r ) + (y 1 1 \...,y 1 Wl \...,y r 1) ,...,y^ Vr) ) (28) 
v v ' v v ' ^ v ' v ' 

JVi N r JVi N r 

Formally, this can be seen as a change of integration variables, where the y-variables describe fluctuations near 
the saddle point. Doing this change of integration variables, we obtain the following integral: 

z = fl fj / d y« ] ex p {-\ v ^ + y^)) x 

x ii n iv^-y^T- n nn(^-^+^ ) -^ ) ) 2 " (29) 

a=l l<i<j<N a l<a<b<r i=l j=l 



G 





Figure 1: Filling numbers Ni describe distribution of eigenvalues between critical points fj,i of V{z). 



Expanding the potential into Taylor series, we find 

* = n n / ^ ] -p f -\n^) - \ £ ^^ {k) 0*0 ) >< 

a=li=l J \ k=2 



r „„ N a N b „„ 

xii n {vp-vp) ■ n nn(^-w+^-^) ( 3o ) 

a=l l<i<j<JV a l<a<6<r i=l j=l 

This integral can be written as 

Z = const x [] II / ^ ex P (4 E %T-^ fc) 0*0 x 
o=li=l"' V fe=2 ' / 



xn n (vf-y^Y- n nn (31) 

a=l l<i<j<A'a l<o<6<ri=lj'=l \ pa p / 

where the normalisation constant is defined as 

const = n exp U*±V{pJ\ ■ J] (f** - ( * b ) 2f,N ' Nl (32) 

a=l ^ ' l<a<fc<r 

The variables . . . can be considered as eigenvalues of iV a x N a matrices $ a , where a takes values 

1, . . . , r. In terms of <J> a , the integral is slightly simplified: introducing the /3-deformed Vandermonde measure 

_ 20 Na r 

Dp*a = LI - ya } ) t rfC/ ] II / d Va 5 . where $ ° = U diag^ . . . /7+ (33) 

and introducing a convenient trace-like notation for the symmetric combinations of variables 

N a 

£(#) fe = Tr(<i> a ) fc (34) 
i=l 
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we can represent the measure (including its potential part) in these new matrix variables: 



r N a ( 1 oo I (i)\k \ r 

nn/^-p -S^w n n (^-^ 

a=li=l J \ k=2 ' / a=l l<i<j<N a 



<3 



2/3 



= f[ / U^exp f-^E^f^^)(, a) ) (35) 

a=l J \ k=2 ' / 

The second factor, which mixes different matrices with each other, can be also represented in terms of traces. 
To do this, one needs to represent this factor as exponent of a logarithm 

N « N * ( ,.w Jj)\ 20 ( N a Nb / (i ) or 

n nn(i+"- -4' s EE'4 + r^ 

l<a<6<ri=lj=l \ Ma Mb / ^ l< a <b<r i=l j=l V Ma Mb 

then formally expand the resulting logaritm into Taylor series 

/ N a N b / „M_„Cj)\\ / JV. JV fc oo f _ nfc (Ai)_..{3)\ k 

y l<a<b<r i=l j = l \ ' ) \ l<a<b<r 1=1 j = l fe=l \ / 

and finally convert the symmetric combinations of variables yi l ^ into traces of matrix variables $ a : 

\ l<„<Kr t ,l,= lfcl K \ Ma ^ b / 

/ E y (-lF( P + g -l)!Tr(d> a ) PT r(^ 

After these algebraic transformations, the partition function is rewritten purely in terms of matrix variables: 
Z = const x J Dp*!... J D P * T cxp £ £ W fe > (Ma)^^ - 



-2/? V V (-iHP + .-^Tr^Tr^, 
It is also convenient to rescale these matrix variables in the following way: 




, . where A a = y( 2 )( Ma ) = y"( M a) (37) 

^a 

Such a rescaling makes the coefficients in front of the Gaussian terms Tr (<J> a ) 2 equal to unity, what is convenient 
for the subsequent perturbative expansion. Ultimately, we obtain an r-matrix Gaussian integral 



Z = Zo x/D^..Y D ^exp(-^-...--|^-^ 1 ,...,^)) 



(38) 



with peculiar interaction 



^ 1 ,...i r )=^^4WTV(i) fc + J2 E B pq (n a ,» b )Tr& a yTr(cf> b y (39) 

a=lfc=3 Ko<Krp+5>0 



where 



h k\ 



A a 



k/2 



B pq (Ha,fJ-b) = 2/3 

and normalisation constant 



i-irip + q-iy. 1 



p\q\ 



h V /2 ( h 



(3Na(N a -1) + N a 



n ^-^?" n ° 



(40) 



(41) 



(42) 



Ka<b<r 



One can see, that at h — > the potential is suppressed by powers of ft. Thus, in accordance with general 
quasiclassical prescriptions, the partition function Z can be treated perturbatively, as a Gaussian integral 



Z = Z x y . . . J Dp4> r ^1 + ^+ ^- + ^- + ...^ cxp 



Tr (0O 2 



Tr {<p r f 



(43) 



Substituting the expression (39) into the expansion (43) and gathering terms of equal ft-degree, we find 




(44) 



This formula is, basically, the main output of the DV method: evaluation of the multiple non-Gaussian integral 
(24) in terms of power series in h. It has three consistutients. The first is normalisation constant Zq, which 
represents the semiclassical contribution (the value of the integrand on the selected saddle point, i.e, on the 
selected solution (Z\, . . . , Zjy) of equations of motion). The second is a product of /3-deformed volumes 



VoL3(A0 = J Drf exp = (V^) N f[ 



NxN 



r(i + ^fc) 



(45) 



of the unitary groups U(N 1 ), . . . , U(N r ). This product is the contribution of the terms of order W° in (43), 
and all the higher contributions are proportional to this product. It is necessary to note, that Vo\p(N) has 
the meaning of the volume of the unitary group only for f3 = 1; for generic /?, we are not aware of any simple 
group-theoretical representation. Finally, the third contituticnt is the quantum series 

Z h = 1 + H k h k 

k>0 

which represent perturbative in h quantum corrections to the leading (nonperturbative in H) classical part Zq. 
Each particular coefficient Hk can be straightforwardly calculated: for example, we find 



1 x . V^{f, a )C 33 (N a ) 1 x . V^{^ a )C,(N a ) x . VW{» a )C 13 (N a )C (N b ) 



72 



E 



24 



E 



A 2 



+ ?E- 



3 ^ A^-^) 



+ 



E 

a#6 



2/3 2 C 11 (N a )C 00 (N b ) + pC 2 (N a )C (N b ) 



2f3 2 C 11 (N a )C (N b )C (N c ) 



A a (^ a - /J, b ) 2 

where Ck lt ...,k m ar e the Gaussian correlators 



A a (n a - 



{fJ-a - Mfe)(Mo - fJ-c) 



/ D^Tr(^...Tr ( ^ exp (-^^) 



Vol /3 (iV) 



NxN 



(46) 



(47) 
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with /3-dcpendcnt eigenvalue measure. We compute these Gaussian correlators, using the Ward identities for 
the integral (47), which can be written in the form of recursive relations 

m 

Cfci,...,fc m = (i ^ Cfei,...,fe m _i,a,6 + ^2 ^ki,...,k t +k m -2,...,k m -i + (1 — P)(k m — l)Cfc 1 ,...,fe m -2 (48) 
a+6=fc,„-2 (=1 

Since each step of induction lowers the total sum of indices k\ + . . . + k m by 2, Ward identities allow to calculate 
each particular correlator Ck lt ...,k m (N) in no more than (fci + . . . + k m )/2 steps. This gives a convenient and 
simple computer algorithm to calculate these correlators. Several simplest correlators for even k\ + . . . + k m are 



level 2 : C 2 = (1 - 0)N + 7V 2 /3, Ci,i = N 



level 4 : C 4 = (3 - 5/3 + 3f3 2 )N + (5/3 - 5(3 2 )N 2 + 2N 3 p 2 , C 3 ,i = (3 - 3/3)7V + 3iV 2 /3 



C 2 , 2 = (2 - 2j3)N + (1 + P 2 )N 2 + (2/3 - 2/3 2 )/V 3 + N 4 (3 2 
C 2 ,i,i = 2N + (1- fi)N 2 + 7V 3 /3, Ci,i,i,i = 37V 2 



level 6 : C 6 = (15 - 32/3 + 32/3 2 - 15/3 3 )/V + (32/3 - 54/3 2 + 32/3 3 )JV 2 + (22/3 2 - 22/3 3 )JV 3 + 5JV 4 /3 3 
C B ,i - (15 - 25/3 + 15/3 2 )/V + (25/3 - 25/3 2 )7V 2 + 107V 3 /3 2 



C 4 , 2 = (12 - 20/3 + 12/3 2 )JV + (3 + 12/3 - 12/3 2 - 3/3 3 )/V 2 + (8/3 - 7/3 2 + 8/3 3 )JV 3 + (7/3 2 - 7/3 3 )/V 4 + 2JV 5 /3 3 

C 3 , 3 - (15 - 27/3 + 15/3 2 )/V + (27/3 - 27/3 2 )/V 2 + l27V 3 /3 2 

C 4 ,i,i = (12 - 12/3)/V + {3 + 7(3 + 3(3 2 )N 2 + (5/3 - 5/3 2 )iV 3 + 2/V 4 /3 2 

C 3 , 2 ,i = (12 - 120)N + (3 + 6/3 + 3/3 2 )/V 2 + (6/3 - 6/3 2 )iV 3 + 37V 4 /3 2 

C 3 ,i,i,i = 67V + (9 - 9/3)iV 2 + 97V 3 /3 

C 2 , 2 ,2 = (8 - 8/3)N + (6 - 4/3 + 6/3 2 )iV 2 + (1 + 9/3 - 9/3 2 - /3 3 )iV 3 + (3/3 + 3/3 3 )/V 4 + (3/? 2 _ 3/3 3 )iV 5 + N 6 f3 3 
C 2 ,2,i,i = 87V + (6 - 6/3)iV 2 + (1 + 4/3 + (3 2 )N 3 + (2/3 - 2/3 2 )iV 4 + N 5 f3 2 
^2,1,1,1,1 = 12iV 2 + (3 - 3/3)/V 3 + 3iV 4 /3, Ci,i,i,i,i,i = 157V 3 



For odd k\ + . . . + k m , correlators vanish due to reflection symmetry of the Gaussian action Tr ((f)) 2 . Note, that 
the above formulas describe correlators with non- vanishing indices ki > 0: every vanishing index fa — simply 
contributes an overall factor of N to the correlator. This (formulas (44) and (42), (45)), (46)) completes our 
description of DV expansion for generic potentials V(z). We now turn to discussion of the particular case of 
3-Penner potential V(z) — a\ log(z) + a 2 log(-J — q) + a 3 log(z — 1). 
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4 /i-expansion for the 3-Penner case 

We are interested in applying the general DV method, described in section 3, to the particular partition function 

"1^2 0^3 Nl q Nl+N2 1 

z DF = q ^ 2 (i-g)wn fdzi n I** n^-^n^ i/ft (^-9r 2/ft (^-ir 3/fi (49) 

i=l q i=N!+l J Q i<j i 

i.c, to the particular log-potential 

V(z) = -on log(z) - a 2 log(z - q) - a 3 log(z - 1) (50) 

This partition function has a DV expansion in powers of ft, which we denote as Z^: 

Z {DV ^ = Z u{1) Z Z h (51) 
Here the first part is a simple additional prefactor from the definition of Dotsenko-Fateev partition function 

Q\OL 2 CX 2 Ct 3 

^ (1) =,W(1-,)W (52) 

while the second (the quasiclassical part Z ) and the third (the quantum part Z h ) are given by formulas from 
section 3 and depend on the critical points - solutions of the equation of motion 

dz z z — q z — 1 

This equation of motion has just two solutions z = ii\,H2'- 



1 (an + a 2 ) + (ai + a 3 )q 1 J{a.\ + a 2 ) 2 - 2{a\ + u\u 2 + aia 3 - a 2 a 3 )q + {a[ + a 3 ) 2 q 2 

Mi,2 = o ; ; T n ; ; v 54 ) 

2 ai + a 2 + a 3 2 a\ + a 2 + a 3 



Accordingly, the quasiclassical part takes form 



PN a (N a -l) + N a 



Z = (Mi - »2) 20N ^ J! VolM^)exp [~V(^ (£) 2 (55) 

and the quantum part takes form 

Z h = 1 + H ^ k (56) 

fc>0 

where the quantities Ai, A2 and Hk are given by formulas from section 3 and straightforwardly computable. 
The quasiclassical part is leading in the regime ft — > and non-perturbative in ft, while the quantum part has 
a form of power series expansion and is suppressed in the limit ft — > 0. However, we are interested not in the 
DV expansion itself, but rather in comparison with a different expansion in powers of q. For this purpose, we 
need to treat all the quantities of interest (including the critical points fi a , the second derivatives A a and higher 
derivatives) as series in q. It is easy to see, that g-expansions for the critical points 



0>1 OL\OL 2 OL 3 2 

/Ji = ■ q - ■ rg q" + ... (57) 



ot\ + a 2 a 2 a 3 a. 1 a. 2 a. 3 2 

H 2 = h 7 77 7<7+7 tt9 +••• ( 5 °) 

«i + a 2 + a 3 (ai + a 2 )(ai + a 2 + a 3 ) (ai + a 2 ) A 
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imply the following expansions for the three parts of partition function: 



Zu(i) = <Z d ° Sl x 1 + J2 E 

\ k—1 m— — k 



2m 



(59) 



oo 



z a = c x h dc ^ x q dc so x i + h ^q k h n 



k—1 m= — k 



(60) 



(61) 



k=0 m=l 

Accordingly, the full partition function posesses the following double (q, h) expansion expansion: 



oo oo 



Z {DV ^ = CX H dcS » X q dcS " X Yl Z kmQ k fi rl 



(62) 



Vfe=0 m=-2k 



Complete information about the partition function is encoded in the following data: the overall normalisation 
constant C, the overall degrees deg^,deg g and the infinite set of coefficients {Zk m }. All these quantities are 
straightforwardly computable with the quasiclassical DV methods of section 3. The constant C is given by 



logC = ^logVol^Ay 



a=l 



— 2^~ Nl + ) ^°g( a i a 2 a 3) + ^Ai(ai logai + a 2 \oga 2 + a 3 log 0:3) 



+ ( 2/3N.N2 + ^N, + %N 2 X + ^—^N 2 + ^A 2 2 -\{cti + a 2 )(N 1 - N 2 ) ) log(ai + a 2 ) 



( 3 — 3/3 3/3 1 

- f 2/3 Ax N 2 + — Y~ N 2 + +jiui+a 2 + a 3 )N 2 ) log(ai + a 2 + a 3 ) 



The overall degrees are given by 



(63) 



59 2/3ft 2 

The first several coefficients Zkm are given by 



A 1 - P AT j_ 1 - P AT P AT* , P A r 2 





fc = 


fc = 1 


m = -4 








m = — 3 








TO = -2 





an 


TO = — 1 





diiCoi + 611 


TO = 


1 


anc 02 + 611C01 + 610 


TO = 1 


coi 


en + &11C02 + anc 3 + fciocoi 


TO = 2 


C02 


C12 + aiic 04 + 611C03 + &ioCo2 



(64) 
(65) 
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^krn 


fc = 2 


to = -4 


^22 


TO = -3 


«22Coi + aii&n 


TO = -2 


O22C02 + ail&llCoi + Glll&io + &22 + »21 


TO = — 1 


ancn + an6nCo2 + &22Co3 + &21 + coian&io + c m b 22 + c ia 2 i 


TO = 


anci2 + &11C11 + a 22 c 4 + au&nCo3 + co2aii&io + C02&22 + c Q2 a 2 \ + &21C01 + &20 


TO = 1 


C21 + &11C12 + hocu + anci3 + a 22 co5 + aii&nCo4 + Q^an^it) + C03&22 + c a3 a 21 + b 2 ic 02 + b 20 c i 


TO = 2 


C22 + &10C12 + »11C14 + 611C13 + &20C02 + &21CQ3 + C04&11&10 + CQ4&22 + C04G121 + G122C06 + ail&llC05 



The first several coefficients a^m are given by 

1 a 2 a 3 1 a 2 a 3 1 a 2 a 3 1 a 2 a 3 1 a 2 a 3 

ail= "2^' ° 21 = _ 4 — ' ° 31 = ~6 — ' a41 = -8— a51= ~10^' (66) 

\a\a\ \oi\a\ VYa\a\ h ol\ol\ 

a 22 =g-^, a 32 =g-^, a42 ---^-, a B2 = (67) 

_ 1 a 3 a 3 1 a 3 a 3 7 a 3 a 3 

° 33 " ~48~ a43 = _ 32^T' a53 = ~192^ (68) 

1 4 4 1 4 4 

1 a 2 a 3 1 a 2 a 3 

fl44= 384^' ° 54 "loO*- (69) 
a55 ~~3840~ (70) 

The first several coefficients bk m are given by 

bio = 7 x (a 2 a 3 7Vi + a 2 iV 2 — aia 3 Ni + aia 2 N 2 — a 2 a 3 Nij3 — a 2 N 2 f3 + aia 3 Nij3 — aia 2 N 2 j3 — 

(ai + a 2 y V 

2a 2 a 3 N 1 N 2 j3 + a 2 a 3 N 2 (3 + ajN^fl + 2a 1 a 3 N 1 N 2 /3 - ona 3 N?P + aia 2 iV 2 2 /3 + 2a 1 a 2 N 1 N 2 /3 + 2a\N 1 N 2 l3^ 

a 2 a 3 N 2 + a\N 2 + axa 3 Ni + axa 2 N 2 

On = ; 

ai + a 2 

x (iaiala 3 Nll3 + a\N 2 + ia\a 2 a 3 N^(3-Aaia 2 alNlj3 2 -\2aia 2 alNiN 2 j3 + 



zu 2( ai +a 2 ) 4 

8a ia2 a 2 3 N 2 N 2 f3 2 - a|iV 2 2 - 4a? a 2 3 N 2 N 2 (3 2 + ha\a\N 2 - 2a\a 3 Ni + 2a\a 2 N 2 + Au x u\N 2 - a\a\Ni - 
a 2 a 2 N 2 - a 2 2 a\Ni - a\N^ 2 - a\N 2 p 2 - a 2 a 2 N 2 - a 2 a 2 N 2 - 2a ia \N 2 - 2a\N^ + 2a 4 iV 3 /3 2 - a^N 2 p + 

a i a 2 a 2 A r 1 2 /3 - aia 2 a|7V 2 /3 + aia 2 a§7V 2 2 /3 - 5aia 2 a 2 iVi/3 + 2a i a 2 a§iV 1 2 /3 2 - 2a\a 2 N 2 (3 - 4a x a^ 2 /3 - 
2aia|iV 2 /3 2 - 4aialiV 2 3 /3 - 2aia 3 ,iV 2 4 /3 2 + 2aia 2 a 2 7V 2 + 5a!a 2 a 2 iVi + a 2 2 a 2 3 N 2 P + 4aiafiV 2 3 /3 2 + 7a 2 a|7V 2 2 /3 + 
2afa 2 iV 2 2 /3 + 8aia^iV 2 2 /3 + 2a\a\Nl0 2 + 3a 2 a 2 a 3 iV 2 + a^a§iVi/3 + cna 2 alN 2 - a\a\Nl(i 2 + a 2 a 2 iVi/3 + 
2af 7ViiV 2 ^ - 2a 2 a 2 7V 3 /3 - a\a 2 3 Nfp 2 + 2a 2 1 a 2 a 3 N 1 N 2 f3 2 + 4aia 2 a 2 , J /V 1 3 ( S - 8aia 2 a 2 7V 1 3 7V 2 /3 2 '+ 
4a^a 3 7Vi^ 3 /3 2 + Saiaaa^TV 2 /? 2 - 4aia 2 i iVi7V 2 ^ - 4aia 2 a 3 7V 1 3 7V 2 /3 2 + 2a\a 2 a 3 N 2 N 2 ft 2 + 2a^a 3 iViiV 2 /3 - 
Aala^NxN^/3 2 - 2a\a 2 a 3 N 1 N 2 2 l3 + a 2 a 2 a 3 A r 1 2 /3 - 2a 2 a 2 a 3 iV 1 2 iV 2 /3 2 + 2a 2 1 a 3 7V 1 2 iV 2 /3 2 + 2a 2 a^iViiV 2 /3 - 
2a|a 3 A^i^V 2 2 /3 2 - 4aia§A^ 1 2 iV 2 /3 2 - 4a? a 3 N 2 N 2 p 2 + 4a|a 3 iVi7V 2 /3 + 2a|al7Vi7V 2 ^ + 4a?a 3 iViiV 2 /3 + 
4a?a 2 A^iA^ 2 /3 - 2afa 3 iViiV 2 /3 2 + 4a?a 2 iViiV 2 2 /3 2 - 4aiaia 3 A^ 1 2 7V 2 ^ - 4a?a 2 7Vi7V 2 3 /3 2 - SafaaTVf^ 2 /? 2 + 
4a 2 a 2 .7V 1 3 iV 2 /3 2 + 2a?a 2 a 3 iViiV2 - 4a?al7V 1 2 A^ 2 2 /3 2 - 4a?a 2 iViiV 2 2 /3 - 4aia§7V 1 2 iV 2 2 /3 2 + 2a 2 a 2 iViiV 2 /3 + 
4aia|a 3 Af 1 2 7V 2 /3 2 + 4aia^iViA^ 2 2 /3 2 + 8a?a|AfiA^ 2 2 ^ 2 - 12a?a 2 a 3 A r i7V 2 /3 - 12aia|a 3 7ViiV 2 /3 - 3aia^a 3 iV 2 /3 + 
3atiO%a 3 NfP - 3aia^a 3 iVi/3 - a?a 2 a 3 iVi/3 + 2a i a 2 a^iV 1 4 /3 2 + Saia^asiV^I/? 2 - 4aia|^Vi^V|^ 2 - 
8aia 2 a§iV 2 7V 2 /3 - 2ala 3 7V 1 2 A^ 2 2 ^ 2 - 4a 2 a 2 N 1 2 iV 2 2 /3 2 - 8a?aliV!iV 2 3 /3 2 - a\a 2 2 N%(i 2 - 2a|a§iV 3 /3 - 
4afN 2 N^/3 2 - 2a\a 2 2 Nlfi - a 2 2 a\N^ 2 - 2a\a 3 N 1 N 2 + Aa\a 3 NlN 2 p 2 + 2a 2 a 2 a 3 N 2 N 2 /3 + 4a?a 3 A^ 1 2 7V 2 /3 - 
8ala 3 N 2 NiP 2 - 2ala 3 N?N 2 /3 + Aa 2 2 a 2 3 NlN 2 p 2 + Aa 2 2 a 2 3 N 2 N 2 l3 + 4a 2 a 2 7V 1 2 iV 2 /3 - 8a 2 a!iViiV 2 2 /3 + 
2a 2 a 2 a 3 7Vi7V 2 /3 2 - 3a?a 2 a 3 7V 2 /3 + 3a 4 iV 2 /3 + 2a?a 2 7V 3 /3 2 + a\a 2 3 N 2 p - a 2 a 2 N 2 (3 2 + 3aia 2 a 3 N 2 + 

a\a 2 a 3 N 1 + 2a 3 1 a 3 N 1 f3 - 2a\a 3 N 2 p + Za 1 a 2 a 3 N 1 + 2a 2 a 2 iV 3 /3 2 - a 2 a 2 N 2 /3 2 - ha 2 a 2 N 2 ^j 
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f 2aia 2 .a 3 iY 2 2 /3 + a\N 2 - 2a 1 a 2 alN 1 N 2 j3 - 2ajN 2 2 + iu\ot 2 2 N 2 + afa 3 Ni + a\a 2 N 2 



Z1 2{ ai +a 2 y 

3atalN 2 + 2a\a 2 3 N\ - 2a\a\~Nl - Aa^a\ iVf - 2a\N%P + a\a 3 N 2 - 2o3 l a 3 N 2 2 - 2aia 2 2 a 3 NiN 2 + 
AalalNxN^P - 4aia 2 a§iViiV 2 2 / 3 + 2a 1 a 2 alN 1 N 2 - 6aia|a 3 A^iA^|/3 - 2aia|a 3 iV : 3 /3 + 2u 1 u 2 u\Nlfi - 
Aa x a\ N$/3 - 2aia 2 alNf - a^ajNx + 2a\a\N%fj + Aa^a^N^ + 2a\a 2 a 3 N 2 + a 1 a 2 a\N 2 + 2a?a§./V 3 /3 - 
2aia 2 a§7V 1 3 /3 - Aa^N^P + Aala^N^ - 10a?a 2 a 3 iViiV| /3 + 2ala 3 N 1 N 2 /3 + 2ala 2 3 N 1 N 2 /3 - 
4a?a 2 ./ViiV 2 2 /3 + 2aia^a 3 7Vi7V 2 /3 + 6aia 2 alNfN 2 P - 2a\a\Nlfi - 2a\a 3 N x N 2 - 2a\a 2 a 3 NlN 2 p - 
Aa\a 3 N 2 N 2 p - 2a z 2 a 3 N 2 N 2 p - 2a 2 2 a 2 3 N 2 N 2 (3 - Aala 2 3 N 2 N 2 f3 - 8a 2 a 2 .7Vi7V 2 £ + 2af TVf /3 - 2a x a\a 3 N 2 - 

2a 2 a 2 N 1 N 2 + 2a\a 3 N 2 fi - 2a 3 a 3 iV 3 /3 - 2a\a 2 3 N 2 (3 + ia ia 2 2 a 3 N 2 + a 2 a 2 a 3 N^j 

(a 2 a 3 N 2 + a\N 2 + a x a 3 Ni + aiu 2 N 2 ) 2 
o 22 = — - 



2(ai + a 2 ) 4 
The first several coefficients Ckm are given by 



coi = — 7 T7 r x a?a?7Vi -3a?a^, 2 /3 2 + 2a2 a 27v3^2 +Q ,3 N Y ~2\a. 1 u\a 3 N 2 $ - 

12a 3 a 2 ai(ai +a 2 )(ai +a 2 + a 3 ) V 1 3 1 3 1H 1 3 1H 2 2 2H 

3a!a 2 alN 2 (3 + 3aia 2 a§.2Vf /3 + 2aia 2 a^iV|^ 2 - 3aia 2 alN$ /3 2 - llaia^iVi^ 2 + 21a 1 a 2 alN 1 f3 + 
2l ai a 2 a 2 N 2 f3 2 + 2\a x a 2 2 a 3 N 2 2 + 2a^a|iV 3 /3 2 - 6a?a 2 ! Ar 2 2 /3 2 + 2a? a 2 ,^ 2 + 2a z 2 a 3 Nlp 2 + 3a 3 a 3 7V 2 £ - 
10aia|a 3 A^i - 3aia 3 7Vi/3 - 3a?a 3 A^ 2 /3 2 - 3a?a 3 iVi/3 + 2a? a 3 NfP 2 - Sa^Nifi 2 + Za 2 2 a 2 3 NtP + 2a x a\Nl0 2 - 
6a? alN 2 p - 3aia|A^ 2 /3 - 3a?a 2 iV|/3 2 - 3a?a 2 N 2 fi + 4alalNl(3 2 + 3ala 3 N?fi + 2a?a 2 iV 2 3 /3 2 - 18a?a 2 a 3 A^ 2 /3 + 

3a?a 2 iV 2 /3 + a 2 a 2 iVi + 24a?a 2 a 3 iViiV 2 ^ 2 - lla?a 2 a 3 iV 2 ^ 2 - Za 2 2 a 2 3 N 2 j3 2 - 3a?a§iVi/3 - 3a|a 3 iV 2 /3 2 + 
a|a|iVi/3 2 + ai alN 2 (i 2 + a\a 2 N 2 l3 2 + a\a 3 N^ 2 - 3a|a§iVi/3 + a 3 a 3 iVi£ 2 + ajajN^ 2 - lla ia la 3 N 2 /3 2 - 
%a ia l<x 3 Nfl3 2 + 21aia 2 .a 3 7V 2 /3 + 24a :L a 2 .a 3 7V :L 7V 2 /3 2 _ lOa^a^/? 2 - 24a?a 2 a 3 iViiV 2 2 /3 2 - 24a?a 2 a 3 /V 1 JV 2 /3 - 

10aia|a 3 iV 2 3 /3 2 + 21a?a 2 a 3 iV 2 2 /3 2 + aia 2 a§iV 2 /3 2 - 8a?a 2 a 3 -/V 3 /3 2 + 18a?a 2 a 3 iV 2 /3 2 + 18a?a 2 a 3 ./Vi/3 - 
10a?a 2 a 3 7Vi^ 2 - lSaia^-a-jiV 2 /? + 18aia|a 3 iVi/3 - lOa^a 2 ^ 3 /? 2 + 24aia 2 a|A^ 2 iV 2 ^ 2 + 12aia 2 a|iViiV 2 /3 - 
21a?a 2 a 3 A^ 2 2 /3 - 10a?a 2 a 3 A^ 2 3 ^ 2 + lSaia^A^ 2 /? 2 - 24aia|a 3 iViiV 2 2 /3 2 + 21a?a 2 a 3 7V 2 /3 - 12aia 2 a|iViiV 2 /3 2 + 
12aia 2 a|7Vi7V|/3 2 + 2a?a|iV 2 + a?a 3 iVi + a\a 2 N 2 + a x alN 2 - 21aia 2 a^ 2 ^ - llaia 2 a 2 .7Vi + 6a?a|A^ 2 2 /3 + 

3aia^|/3 - lla?a 2 a 3 A^ 2 + aia 2 a^iV 2 - 24aia 2 a 3 iViiV 2 /3 + 3a?a^iV 1 2 /3 - llaia 2 a 3 iV 2 - 10a?a 2 a 3 iVi^ 

c n = - 2Nl( > ai ~" 2) x fa 3 - 2a 3 ^ + a 3 f3 2 - a 3 N 2 (3 + 2a 3 N^ + a 3 N 2 /3 2 - 2a 3 N 1 /3 2 - Za^Nrf 2 + 
a 3 N 2 p 2 + a 2 N 2 p 2 - a^N^ 2 + ai N 2 p 2 ~ a^Nrf 2 } 

c 2 i = — x f 12aia|iViiV 2 2 ^ 2 + 20a?a^iVi7V 2 2 ^ 2 - 6a^a?7ViA^ 2 /3 2 - 6a?a§iVi + \2a\a\Nl0 2 - 

2(ai + a 2 )° \ 

6a?a§7V 1 3 ^ 2 - 2a 3 ,a 3 A r i - haiala 3 N^P + aia 2 alN 2 j3 - a 1 a 2 a§iV| J 8 + aia 2 a§A^|/3 2 + 21aia 2 a§7Vi^ 2 - 
41aia 2 a|7Vi/3 - 41aia 2 a|7V 2 /3 2 + 5aiala 3 7V|/3 2 - ealaiiVf/3 2 + 8a?a!iV 2 2 /3 2 - 4a?a|iY 2 /3 2 - 2a|a 3 A^ 1 3 /3 2 - 
12a?a 3 ^V 2 /3 + 9aia|a 3 iVi +4afa 3 iVi/3+ 12a?a 3 A^ 1 2 /3 2 + 12a?a 3 A^i/3- 6a?a 3 iV 3 /3 2 +4aiai7V 2 2 /3 2 - \2a 2 2 a 2 3 N 2 [3 - 
2a 1 a\Nlfi L + 8a?a|7V 2 ^ + 4aia|iV 2 /3 + 4a?a 2 iV 2 2 /3 2 + 4a?a 2 7V 2 ^ - 4a?a|iV|/3 2 - 4a|a 3 7V 2 ^ - 2a?a 2 iV|/3 2 + 

9a?a 2 a 3 7V 1 2 /3 - 4a?a 2 AT|/3 - 6aia§A^i + 2a?a 2 a 3 iVi7V 2 /3 2 - 3a?a 2 a 3 iV 2 /3 2 + 12a|a|iV 2 ^ 2 + 12a?a§7Vi/3 + 
4a|a 3 7V 2 /3 2 - QalalN^ 2 - 2a x alN 2 p 2 - 2a\a 2 N 2 p 2 - 2a 3 a 3 7Vi/3 2 + 12a^a|iVi/3 - 6a?a 3 7Vi/3 2 - 6a?a^7Vi/3 2 - 
3aia|a 3 A^ 2 /3 2 + 8aia^a 3 iV 3 /3 2 + 5aia^a 3 7V 2 /3 + 6aia|a 3 iViiV 2 /3 2 + 9aia|a 3 7Vi/3 2 + 30a?a 2 a 3 7Vi7V 2 2 /3 2 _ 

2a?a 2 a 3 iViiV 2 /3 - 2aia|a 3 A^"|/3 2 + 5aja 2 a 3 N^l3 2 - a x a 2 alN 2 l3 2 + 4a?a 2 a 3 iV 3 /3 2 - 9a?a 2 a 3 iV 1 2 /3 2 - 
9a?a 2 a 3 iV 1 /3 + ha 2 1 a 2 a 3 N 1 l3 2 + 17a 1 a|a 3 7V 1 2 /3 - 17aia^a 3 7Vi/3 + 20aia 2 a^ 3 ^ 2 - 74a!a 2 a^iV 2 iV 2 ^ 2 - 
20010201^1^2/3 - 5a?a 2 a 3 iV 2 2 /3 - 2afa 2 a 3 N^(3 2 - Yla. x u 2 2 a 3 Nl{3 2 + 30aia|a 3 AfiAf|^ 2 + 5a?a 2 a 3 iV 2 /3 + 
20aia 2 a|A^iAf 2 /3 2 + Ma^a! A^/? 2 ~ 4a?a 2 7V 2 - 6a?a 3 iVi - 2a?a 2 7V 2 - 2aia^7V 2 + Ala x a 2 a 2 3 Nl ft + 
21aia 2 a^7Vi - 8a?a|7V|/3 - 4aia 3 ,7Vf/3 - 3a?a 2 a 3 7V2 - aia 2 a^7V 2 + 2a^a 3 7Vi7V 2 /3 + 6a|a§7Vi7V 2 /3 - 
6aia^a 3 7Vi7V 2 /3 - 4afiVi7V 2 2 /3 2 + 4a^ 2 7V 2/ 3 2 + 20a?a§A^i 2 iV 2 ^ 2 + 2aia^iViiV2/3 - 4a?a|iV 1 iV 2 ^ 2 - 
2aia|iViiV 2 2 /3 2 - 6a|a 3 AfiA^ 2 2 /3 2 + 20a 2 a|^i 2 7V 2 ^ 2 + 12aia 3 7V 2 iV2/3 2 - 6a?a|iViiV2^ 2 + 4a?a 2 iY!iV2/3 - 
36a?a 2 a 3 7Vi 2 7V 2/ 3 2 - 8a 2 alN?N 2 j3 2 - 48aia|a 3 iV 1 2 7V 2 /3 2 + 24a?a 3 A^i 2 iV 2 /3 2 + 6a?a 3 iVi/V 2 /3 - 6a?a 3 iViiV 2 ^ 2 + 

2a?a 2 iViiV 2 /3 - 2a|a 3 A^i^ 2 /3 2 - eafagiVxiV 2 /? 2 - 2a?a 2 iV 1 /V 2 /3 2 + 2a\a 2 NlN 2 (3 2 + 4a?a 2 A^iAT|/3 2 + 
6a?a|7Vi7V 2 /3 - 12a?a§iV 2 ^ - 3aia|a 3 Ar 2 + 5a 2 a 2 a 3 /Vi - 2aia 3 /ViiV2/3 2 - 6aia 3 ^ 2 7V 2 /3 2 - 2a 2 a|7Vi7V 2 2 ^ 2 



These somewhat lengthy expressions describe the first several orders of double (ft, q) expansion of Zuv- We now 
proceed to derive this expansion from a different point of view: namely, starting from the contour-integration 
method and the corresponding conformal-ficld-theory g-expansion. 
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5 g-expansion for the 3-Penner case. Comparison 

On the other side, the partition function (49) is known [5, 6] to posess an expansion in powers of q 

z ic FT) = qdcg x ^ ^ («i+a, + 2/3iVi5 f ) x (l + g B fc( /) (71) 

with the overall degree 



fe=i 



<l<^ = ^+^1(^1-1)+^ + ^^^ ,72) 



the normalisation constants 

_ r(x + 1 + /3(k - i))T( y + 1 + p(k - i))r(i + gfe) 

CJV(X ' V) 11 r(.x + y + 2 + (A r + fc - 2)/3)r(/3 + 1) (7d) 

and the first coefficient 

Bi = - (2/3Ar a - 2iVi/3 2 + 2/3 2 J /V 1 2 + 2JVi/3ai + 2a 2 - 2f3a 2 + 2j3N ia2 + a 2 a x + a\) x 

x (4AWi/? 2 + 2/3N ia3 + 2(3 2 Nl - 2(3 2 N 2 + 2f3N 2 a 2 + 2f3N 2Ctl + 2N 2 a 3 fi + 2(3N 2 + a 3 a t + a 3 a 2 ) x 

x (2f3y 1 (2(3N 1 + ai + a 2 )- 1 (2(3N 1 + ai + a 2 - 2/3 + 2)" 1 (74) 

The expansion Z^ CFT ^ was obtained in [5, 6] via exact non-Gaussian contour integration over the contours [0, q] 
and [0, 1] with multiplicities N\ and N 2 , just as indicated in (49). This expansion coincides with standard series 
for the 4-point spherical conformal block, thus we call it a "CFT" expansion. To compare 

Z (CFT) with Z (DV) 

of the previous section, we need to compute the asymptotical ft-expansion both for the Gamma-functions in the 
normalisation constants Cn(x, y) and for the rational expression for B\. This computation is a straightforward 
exercise in undergraduate calculus and reduces essentially to use the Stirling formula 

log T{x) = x log x - x - i log x + 1 log(27r) + ±- - — + O (J^) (75) 

at large x. With dimensional analysis one can easily estimate that Z^ CFT ^ has precisely the same form as Z^ DV ^ 

(GO OO \ 

J2 Z km q k n m \ (76) 
k=0m=-2k J 

with some degrees deg s ,deg g and coefficients C, Zk m , which need to be determined with the Stirling formula. 
Actual symbolic computation with the help of MAPLE shows, that C — C, dcg h = deg ft , deg g = deg g and 

Z km = Z km , at least for (fc, m) - (0, 0), (0, 1), (0, 2), (1, -2), (1, -1), (1, 0) 
thus suggesting that 

~ ! (77) 



Z (CFT) = Z (DV) 



This completes our check of the relation between two different methods to treat the 3-Penner ensemble and the 
two corresponding expansions. The detailed intermediate calculations with the Stirling formula, which lead to 
equalities Z km = Z km are too lengthy to be presented here. We only derive now the first three equalities, of 
normalisation constants C = C and of overall degrees deg ri = deg fi , deg g = deg g . 
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First of all, the overall degrees in q coincide trivially: deg q — deg g = deg, so we are left only with the 
check of coincidence of overall degrees in ft and of the normalisation constants. The small-ft asymptotic Stirling 
expansion, necessary for such a check, includes terms of order ft -1 , ft and log ft: 



logCV, ^ — , —J + logCV 2 



^±^ + 2PN u ^\ =^logVol (Ar a ) + 

' a=l 



log r + 1 + /3(fc - i)) + io g r + 1 + - i)) - log r + 2 + p(m + k-2) 



+ 



N 2 r y 

+ Xj logr^ " 1 ^" 2 +2/3jv 1 + i + /3(fc-i)j + i og r(^ + i + /?(fc-i)) 



_ logT ( Ql + " 2 + 03 + 2/3iVi + 2 + /3(iV 2 + k - 2) 



^logVol^A^) + ( 1 - - A^ + -iVf J log(aio;2a3) + ^./Vi(ai logai + a 2 loga 2 +a 3 loga 3 ) + 

a=l ^ 



- ( 2(3N 1 N 2 + ^-^N 1 + ^-Ni + ^^N 2 + ^Nl-^(a 1 +a 2 )(N 1 -N 2 ) ) lo S (n + n_>) - 



3 — 3/3 3/3 1 
- ( 2^iViiV 2 + — ^"^2 + + - (ai + a 2 + a 3 )N 2 ) log(ai + a 2 + a 3 ) + 



^T^i + + ?JV? + £jV| ) log ft + OCft 1 ) = log C + (deg h ) ■ log ft + Otft 1 ) 



(78) 



i.e, we indeed have deg^ = deg^, deg g = deg 9 , as expected. Similarly, inclusion of higher order terms (ft 1 , ft 2 
and so on) in the Stirling expansion reproduces the coefficients Z i,Z 02 and so on: 



logCW, ^) C N2 + 2/37Vl ' t) = lo & C +( dc &n) ■ logft + Z 01 ft+ (z a2 -3^h 2 + ... (79) 



To reproduce the coefficients Z^ m with k > 0, it is essential to include the coefficients S/j into consideration. 
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6 Conclusion 

In this paper we demonstrated the consistency of two different treatments of the multi-Penner /3-ensemblc, 
which plays an important role in the study of AGT-relations. 

• One side of equality is the original Dotsenko-Fateev integral of [5, 6] with a system of open integration 
contours taken with multiplicities (Ni, . . . , N r ). In the case of r = 2 the answer was found in [5] and summarized 
in s.5 of the present paper, where it is called Z^ CFT \ 

• Another side of equality is the same integral, but treated differently: as Givental decomposition [13, 15] 
into a system of Gaussian models with the filling numbers (Ni, ...,N r ). Above sections 3 and 4 are devoted to 
quasiclassical evaluation of this decomposed integral, in the form of a series in parameters eti, . . . , a r+ \. The 
result for r = 2 is exposed in s.4, where it is called Z^ DV \ 

On the two sides of the equality the integrands are the same, but integration contours are different. Not- 
surprisingly, expansion coefficients in the two cases look quite different, but in fact they coincide(l), as stated in 
eq. (77) . A toy example of this phenomenon - of contour independence of quasiclassical expansion - is analyzed 
in s.2. Significance of the equality (77) between two different expansions 

y{CFT) _ 7 {DV) 



is that to the r.h.s. one can apply the well developed matrix-model technique. In particular, F = log Z^ DV ^ 
should have a representation, typical for the Seiberg-Witten (SW) theory [11]: 

where the role of the Liouville-SW differential is played by exact resolvent 

where the averaging is that of the multi-Penner /3-ensemble. After this, the AGT relation implies that such 
generalized Seiberg-Witten prepotential is logarithm of the ei^-regularized partition function of ADHM super- 
instantons. This identification still remains to be understood. 

In fact, (77) is not fully proved in this paper, only the first terms of the q, ^-expansion are explicitly evaluated 
and compared. However, already this check is somewhat involved. 

The difficult part is Z<- DV K Notably, F = logZ( w ) is a well-known quantity. It is actually the cele- 
brated CIV-DV prepotential [13], just modified (generalized) in three different directions: from polynomial to 
logarithmic potential V, to ft ^ 1 and to arbitrary genus. 

In fact, the genus expansion of matrix model free energy is not the quasiclassical expansion of the sections 
3 and 4. Quasiclassical expansion is in h at given values of N a , while genus expansion is in h at given values 
of S a = HN a . These two expansions are of course related, but different. We devote a big appendix below to 
demonstrate that the same general DV-phase formulas (44), (42), (45)) and (46)), which are used to reproduce 
Z(CFT) £ or logarithmic potential V in s.4, in another case - for the polynomial potential V - reproduce the 
well-known CIV-DV prepotential with above-mentioned modifications. This brings the well-developed theory 
of CIV-DV prepotential into the circle of subjects, embraced by the AGT relations. 
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7 Appendix, /i-expansion for the polynomial case. Generalized 
CIV-DV prepotential 

This appendix is devoted to another special case of general constructions of section 3 - the special case of 
polynomial potentials V(z) of arbitrary degree deg V(z) = r: 



V(z)=T +T lZ + ... + T r+1 z r+1 

As in the generic case, these formulas are conveniently expressed not through the coefficients T a , but rather 
through the critical points fi a - roots of the derivative 

V'(z) = {z-ni)...{z-nr) (81) 

However, the polynomial case is certainly much simpler than the generic case, since all the quantities of interest 
- including the values of potential and its derivatives at the critical points - are completely determined here by 
the critical points theirselves. For example, for the 2nd derivatives we have 

A a = V"(fl a ) = H(»a-»b) 

and expressions for higher derivatives are equally straightforward. These expressions need to be substituted 
into the generic quasiclassical formula for the partition function 

Z = Z Q xf[Vo\p(N a )x ll + J2 H k?i k ) (82) 

o=l \ fe>0 / 

or directly into the free energy 

r f H 2 \ ( H 3 \ 

F = \ogZ = logZo + ^logVoLj(Ag + H 1 h+ {H 2 - -fj h 2 + \ H 3 - H X H 2 + -±J h 3 + ... (83) 

After the substitution, the answer becomes expressed entirely in terms of independent variables of the polynomial 
model: namely, the critical points . . . ,/z r , the filling numbers Ni,...,N r and the parameter (3. However, 
as explained in the Introduction, instead of the variables N a it is conventional to introduce different variables 
S a by the rule S a — hN a . As follows from this rule, the number of integrations N a is no longer a constant: it 
becomes proportional to hr 1 . Quasiclassical expansion with such an additional requirement is usually called a 
genus expansion. Throughout this section we express the free energy in terms of variables S a , thus dealing with 
a genus expansion, not just with quasiclassical expansion. 

As it is easy to find out, the genus expansion of the free energy has a form 



F = /Z [G k/2 (Si,...,S r ) + F k/2 (S 1 ,...,S r )]h k - 2 (84) 



k=0 



with two types of contributions: the so-called perturbative prepotentials Gh/ 2 (Si, . . . , S r ) (they receive con- 
tributions from the normalisation constant Z and the /3-deformed group volume factor) and non-perturbative 
prepotentials F k / 2 (Si, . . . , S r ) (coming from the quantum part of the partition function, i.e, from the coeffi- 
cients Hh). The lowest of these quantities - Go (Si, . . . , S r ) and Fo(S\, . . . , S r ) - are known as perturbative and 
non-perturbative CIV-DV prepotentials, respectively. 

7.1 The case of /3 — 1 and genus zero: ordinary CIV-DV prepotential 

Using the above general formulas, one can easily find the perturbative 3 part of the CIV-DV prepotential, by 
putting /3 = 1 and extracting the contributions of order hr 2 (in terms of variables S a ) from the logarithms of 
both the normalisation constant and the group volume factor. The answer is 



a=l ^ \ a / / a / b 



(85) 



3 Note, that the terms perturbative and non-pcrturbative refer to 4d gauge theory interpretation of these quantities, which does 
not coincide with their matrix-model meaning: in the matrix model, quantities F k ^ 2 (Si, . . . , S r ) come from the perturbative (in the 
small parameter h) part Zn of the partition function, while G k /2(S\, . . . , S r ) originate from the leading non-perturbative (in h) part 
Zq. This interesting phenomenon (perturbative gauge theory quantities correspond to non-pcrturbative matrix-model quantities 
and vice versa) can be considered as a concrete example of gauge-string duality. 
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Similarly, the non-perturbative part of the CIV-DV prepotential can be found, extracting the contributions of 
order hr 2 (in terms of variables S a ) from the logarithm of the quantum part Z%. The answer is an infinite series 
in S 1 - variables. It is convenient to denote its homogeneous parts of degree p as Fq P \Si, . . . , S r ): 



F (S U . ..,S r ) = F^(S U . ..,S r ) + F w (Si, . . . , S r ) + 0(S b ) 



(86) 



For the cubic contribution, we find 



V 



11 15 

a^ c ~ - Mc) Kl 82 ^ + 2 SaSbSc + 2 S « SC + 12 S " 



(87) 



For the quartic contribution, we find 



^ - -3 i 21 2 C 2 | 47 3 



" J^ c A2( Ma -^)3(^-^) ( 2SaS * Sc + l SaS > + 7S " SbSc + l s - Sb + AS " Sc + Y s3aSb + T s « ) + 



A2( Aia - Mb )2( Ma _ Aic )2 (^„S 6 S? + |s„S 6 2 S c + 5^ c 2 + -^S a 2 S b S c + S 2 a S 2 b + 45 a 3 5 c + 



^2 \ ^ \ 



S a SbS c Sd + —S a S b Sd + —S a S b S c 



7 T 5 11 29 

+2S 2 a S c S d + -S 2 S b S d + -SlS b S c + -S^S 2 + 3SlS d + 3SlS c + —S^St, + 



+ £ 



a 9 / T7 T7 77 7 [ Tn'-'aScSdSe + -TzS a S b SdS e + —S a S b S c S e 

Ai(fi a - (i b )((i a - ii c )(p a - ii d )(fi a - fi e ) \12 12 12 



1 5 5 5 5 5 5 

— S a S b S c S d + —S a S d S e + —S a S c S e + —S a S c S d + -^S a S b S e + —S a S b S d + —S a S b S c + 



5 a 5 i ,(4S2 + ii5 a 5 6 + 452) ^ S a S b (S a + S c )(S b + S d ) ^ S a S b (S a + 2S b )(S b + S c ) 



V 



C7 , a 2 Oa " P-bYiHa - Hc)(Hb ~ Md) ^ (Ma - Mb) 3 (M6 - Mc) 



/ 



(88) 



These formulas reproduce and generalise various explicit results, available in the literature [13], in particular 
the cubic contribution coincides with the one obtained in [14] by different methods. 
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7.2 Generalisation to f3 — 1 and higher genera 



To generalise the CIV-DV prepotential to higher genera (i.e, to higher k > 0) one just needs to extract the 
contributions of higher orders h k ~ 2 (in terms of variables S a ). Doing so with the normalisation constant and 
the group volume factor (actually, for j3 = 1 and k > only the group volume factor contributes), one obtains 



and, for all to > 1 



G 1/2 {Si, . . . , S r ) — S a 



o=l 



l()jv(277) - I ( y 



S a 



Gi(s u ..., s r ) = ^ lo s ( y ) + C'(-i) + i lo s( 27r ) 



a=l 



C m +i/2(S'i, • ■ • j S r ) 



G m +i(Si, . . . , 5 r ) 



B 



2m+2 



i=i (2to+ 1)(2to + 2) ^ VS'a 



2m-l 



2m+4 



/3=i (2to + 2)(2to + 4) ^ V5„ 



(89) 
(90) 

(91) 
(92) 



where B/~ are the Bernoulli numbers: 



Bi = 7, B 4 
6 



11 °° 
— , -Bg = — - , . . . , with generating function 



fc=2 



B k z k 
k\ 



— - 1+ 2 



and C(x) = X^i^i * x is the Riemann zeta function. Doing so for the quantum part of the partition function, 
one obtains non-perturbative prepotentials as infinite series in S*- variables: 



F k/2 (S U . . . , S r ) = Ft +3 ~ k \S u ...,S r ) 



(93) 



where Fjf\Si, . . . , S r ) denote homogeneous contributions of degree p. Notably, for j3 = 1 all the contributions 
of half-integer genus (corresponding to odd values of k = 1, 3, 5, . . .) vanish: as we will see in the forthcoming 
section, they are proportional to /3 — 1. All the integer genera do contribute. For genus one, we obtain 



F^Su ... , S r ) = F^(S U . ..,S r ) + F™(S U ...,S r ) + 0(S 3 ) 



(2), 



(94) 



where 



F^\Si, ...,S r ] 



0=1 



£ 



S a 



+ £ 



S a 



±h 6A a (p a - Vb) 2 ' 24A a (yU Q - Hb)(Pa - Mc) 



(95) 



p(2)/Q Q \ _ 1 \ - SgS 



^ajj I \ ~" 1 / 25 7 2 

i(/ia - M & ) 4 ^ A^a - ^) 4 1 12 ° 6 3^ 



+ £ 



a^b^c 



V 



+ E 



20 



+ -TT7 T7 T7 T7 \ { TH^aSe + —S a Sd + —S a Sc + —S a Sb + -TT-zSa) (96) 

H ^ 6 ^ d ^ e A a(Ma -Mb)(Ma ~Mc)(Ma ~Md)(Ma ~Me) V24 24 24 24 288 / 

To obtain explicit formulas for genus two and higher, it is necessary to calculate higher coefficients Hk of the 
quantum partition function Z^. The above results are obtained, making use of coefficients H\ and H2 only. 
Higher Hk can be straightforwardly calculated with the methods of section 3. 

7.3 Generalisation to arbitrary /3 and genus zero 

To generalise the CIV-DV prepotential to arbitrary values of 6, one just needs to keep B as a free unconstrained 
parameter in all the formulas. This slightly boosts the complexity of intermediate algebraic calculations, but 
does not conceptually affect neither the method of calculation, nor the structure of the answer. Extracting the 
contributions of order hr 2 from the logarithms of both the normalisation constant and the group volume factor, 
we find the /3-deformed perturbative CIV-DV prepotential 

'PS 2 (S a \ 3/3, 

a=l ^ \ a/ , a ^ b 



G (Si,...,5 r ) = E( VMaN "T^" WO) -;8£ S « S *lDgk,-w| (97) 



I.e., at the spherical level the /3-deformed perturbative prepotential depends on /3 just linearly. Similarly, the 
non-perturbative part of the CIV-DV prepotential is found extracting the contributions of order fir 2 from the 
logarithm of the quantum part Z%. Direct calculation shows that 

F^{S U . . . , S r ) = P 2 ■ F< 3 \s u . . . , S r )\ (98) 

F (4) (5!, . . . , S r ) = /3 3 • F (4) (5!, . . . , S r )\ (99) 

It is clearly seen, that the perturbative /3-deformed CIV-DV prepotential depends on (3 homogeneously: 

F (AS 1; . . . , XS r ; 0) - AF (Si, ...,S r ; A/3) (100) 
Equivalently, this statement can be expressed in a form of an exact differential equation: 

*f? + *=S s -S (WD 

As we will see, this property fails at higher genera, but possibly gets substituted by some more sophisticated 
and yet unidentified property, which relates the dependence on /3 and the dependence on S'-variables. 

7.4 Generalisation to arbitrary (3 and higher genera 

As already mentioned above, to handle the higher genera (i.e, higher k > 0) one just needs to extract the 
contributions of higher orders H k ~ 2 in terms of variables S a - For generic values of /3, contributions of odd values 
of k (non-integer genera fc/2) become non- vanishing. The simplest of those is contribution of genus 1/2. 

Genus 1/2. For the perturbative prepotential of genus 1/2 we get 

'1, fS a \ B. fS a A a \ l + log(/3)+/3-/31og(/3) 



For the non-perturbative prepotential of genus 1/2, we obtain 



(102) 



F 1/2 (Si, . ..,S r ) = F${Si, ...,S r ) + F${Si, ...,S r ) + 0(S 4 ) (103) 
where the quadratic contribution has a form 
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008-1) f3 „ 3 c2 



E 



(104) 



and the cubic contribution has a form 



-F , j 2 (.s'i,...,,s' r ) — - 



/3 2 (/3-l) /21 2 71 2 59 3 



E 



A2( M a-^) 3 ja-Mc) ( 7 ^^ 5C + 2^ + ^ + i S « Sh + 



- E 



a^b^c 



AK^-l^ydl-^y { SaS2c + T SaSbSc + SaS2b + T 5 ^ c + T 5 ^" + Y s « 



y j A 9 / \o7 \7 \ f %S a S c Sd + —S a SbSd + —S a SbS c + —S a Sfr H TT^a^d + 

a ^ #d A a(Ma - Mfe) 2 (Ma - Mc)(Ma - Md) V 2 2 4 2 



13-. 47 2 49 3 \_ v W-l) 

2 a C 4 «* b+ 4 A^^-^)^-^)^-^)^-^) V24 



~z~;S a SdS e 



17„,„ 17 ~ . 17 „ 17, 



197 



24 



SaS c Se + —SaScSd + —S a Sl,Se + —SaStSd + —SaS^Sc + —S a S e + —S a Sd + — 5 a 5 c + — 5 a S(, + —7 ~5 a I — 



24 



24 



24 



24 



24 24 24 24 144 



E 



/3 2 (/3-l) 

A a A 6 



135 a 5b(5 a + 5;,) SaSt(Sa + Sf, + S c + Sd) 5 a 5;, (25 a + 35j + 25 c ) 



(105) 



Obviously, a variant of the homogeneity property holds here as well: 

F 1/2 (XSi , . . . , AS,. ; /3) = ^§^1/2 (-Si, ... , 5 r ; A/3) 



(106) 



Genus 1. For the perturbative prepotential of genus 1, we obtain 
1 + 3/3 + /3 2 ^,__fS a \ , log(27r) /3 7 



Gi(S'i, . . . , S r ) 



12/3 



a=l 



4 12 12/3 



«'("!) + |||^ (107) 



where 7 = 0.57721 ... is the Euler constant (this is an asymptotic expansion only, its reformulation in terms of 
convergent series or integrals remains to be found). For the non-perturbative prepotential of genus 1, we obtain 



F 1 (S 1 , ... , S r ) = fi (1) (Si, . . . ,S r ) + FF>(S U ...,S r ) + 0(S d ) 



(2), 



(108) 



where the linear contribution has a form 
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fP(s 1 ,...,s t ) = 



Sa 



a^b 



aifJ-a - Mfc) 2 V 6 2 6 



Sa 



a^b^c 

and the quadratic contribution has a form 



11 V±V 



A'&XMa-Mc) V24 8" 24 



a^b 



A2(/x a -/i fe ) 4 



+ 



+ E 



a^b^c 



.115 . 71 H5 ^ \ 

+ (— + — 13 ) 5 ' 



+ E 



a^b^c 



16 ^ 16 ^ + 16 ^ 1 ba 



+ E 

a^b^Cy^d 



AaiVa - Mb) 2 (^a - Mc)(Ma - Md) 



+ 



+ E 



&l(Pa - Mfc)(Ma - Mc)(Mo - Md)(Ma - Me) 



+ 



y— 

^ A., 



a^b 



AaAfc 



- W ~ l) 2 )SgSb 
4(Ma - Mb) 4 



E 



p(p - l) 2 S a S b 



2g(g - l) 2 SgS b 
, a 2(^a - Mb) 2 (Ma - Mc)(A*b - Md) ^ (Ma ~ Mb) 3 (^b - He) 
d=£b 



E 



7 

(109) 



In this case, the prepotential does not seem to enjoy any homogeneity properties. 
Genus 3/2 and higher. For the perturbative prepotential of genus 3/2, we obtain 

G 3/2 (S 1 ,...,S r ) = l + ^L 

For the non-perturbative prepotential of genus 3/2, we obtain 



(HO) 



F 3/2 (Si, . . . , S r ) = (Si, ...,S r ) + 0(S 2 ) 



(HI) 



where the linear contribution has a form 
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2 6^ 6 M 2 P J a ££ c Al(ji a -vL b )3(» a -» c ) 



I , 165 „ 165 ^, , ~\ S a 

16 16 / ^ c A «(M«-M6) 2 (Pa-Mc) 2 



/ 21 111 111 2 21 3 \ ^ S a 



V 16 288 M 288 M 16^ / ^f-J 

' n —f- h —f- f —f- 



To obtain explicit formulas for non-perturbative prepotentials of genus two and higher, it is necessary to calculate 
higher (k > 2) coefficients H^. Much simpler are perturbative prepotentials - they are obtained by expansion 
of the group volume factor and, hence, do not require knowledge of complicated expressions for Hk at all. For 
the perturbative prepotentials of half-integer genera we find 

G ro+1/2 (Si, . . . , s r ) = Q + {2m + i)(2m + 2) § (i) (112) 

For the perturbative prepotentials of integer genera we find 

<W*.-.*> - - ('ffe..- 2 .a. r(28 + 1) r r g: ) _ 28 + 3) g (^) 2 "' (113) 

This completes our description of /3-deformed higher-genera CIV-DV prepotentials. It is important to note, 
that all the non-perturbative prepotentials described above, posess one and the same structure: 

F (3-2g) = y^ H {1A) (S a ,S b ,S c ) H {2) (S a ,S b ) 

9 aib^c Aa ^ a ~ ~ He) f£ Aa(Ma - fib) 2 

p(i-2g) _ ^(1,1,1,1) (S a , S b , S c , Sd, S e ) y-^ H(2A,l)(S a , $b, S c , Sd) 



H( 3S )(S a ,Sb, S c ) ^ H( 2 .2)(Sg, Sb, S c ) \ - H^(S a ,Sb) 

^ Al(p a Mb )3( Mo - /ic ) + Al(n a - ^{Ha He) 2 + ^ A^a " 



a^tb^tc a^tb^tc " a=£b 



y-^ y-^ H(l,l),(l,l)(S a ,Sb,S c ,Sd) 

h> <M*b ^Abifia-fibfiPa-Hc^b-Hd) 



y^ y^ H(2),(i,i)(S a ,Sb,Sc) y^ H/ 2 ),(2)(S a , S b ) , j _ 

f^ b ^ b A a A b (fi a - fi b ) 3 (p b - He) A a A b (fi a - fi b ) 4 
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where #n,i),#( 2 ) are certain polynomials of degrees 3 — 2g in S-variables, while #(1,1,1,1), #(2,1,1); #(3,1): 
#(2,2): #(4)i #(1.1), (1.1) j #(2). (1,1) j #(2), (2) are certain polynomials of degrees 4 — 2g in S- variables. Already the 
examples described above (for genera g = 0, 1/2, 1, 3/2) provide enough evidence to conjecture, that different 
structures appearing in the generalized CIV-DV prepotentials are labeled by (collections of) Young diagrams. 
It would be interesting to check this conjecture and understand, what kind of structures actually appear in the 
next (third) order of S'-expansion. Such understanding could be essential for future applications of the CIV-DV 
prepotentials as yet another special functions of string theory [15]. 
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